ON APPLICATION OF MAXIMUM LIKELIHOOD 
METHOD AND KALMAN-FILTER TECHNIQUE 
TO ESTIMATE PARAMETERS FROM FLIGHT 
DATA OF ROCKETS AND SHELLS 


By 

K. Venkat Chandra Sekhar 


1 % 

cMh 



DEPARTMENT OF AEROSPACE EMOINEERING 


indain Institute of Technology Kapnur 


JUNE, 2004 



ON APPLICATION OF MAXIMUM LIKELIHOOD 
METHOD AND KALMAN-FILTER TECHNIQUE TO 
ESTIMATE PARAMETERS FROM FLIGHT DATA OF 
ROCKETS AND SHELLS 


A Thesis Submitted 

in Partial Fulfillment of the Requirements 
for the Degree of 


Master of Technology 


by 

K. Venkat Chandra Sekhar 



Department of Aerospace Engineering 

Indian Institute of Technology, Kanpur 


June 2004 






2^o 




CERTIFICATE 

It is certified that the work contained in this thesis entitled, "On application of 
Maximum Likelihood method and Kalman-Filter technique to estimate parameters 
from flight data of rockets and shells” by K. Venkat Chandra Sekhar been carried out 
under my supervision and that this work has not been submitted elsewhere for a degree. 




Dr. A.'K/ Ghosh 

Associate Professor 


June 2004. 


Department of Aerospace Engineering, 
Indian Institute of Technology, 

Kanpur -208016. 



ABSTRACT 


An attempt is made to estimate aerodynamic parameters using radar tracked 
coordinates of a free flight artillery shell. The present thesis investigates applicability of 
EBCF and ML method to estimate aerodynamic parameters using flight data with limited 
information content. Due to non-availability of flight data, simulated flight data using 
Point Mass and Six-degree of fieedom model have been generated. 

In particxilar, study was conducted to investigate the applicability of these methods 
in extracting aerodynamic parameters by processing Radar tracked data obtained through 
different class of flight vehicl^store. These methods have been applied to fli^t data 
starting from a one-dimensional ballistic target to, an artillery shell, longitudinal motion 
of Short Range Strategic Missile (SM) and also to flight data of an aircraft. 

The applicability of the existing methods (ML and EKF) is tested for these three 
different classes of flight data. It is observed that both ML and EKF can be 
advantageously applied to estimate few aerodynamic parameters with acceptable level of 
accuracy. 
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CHAPTER 1 

INTRODUCTION 


Artillery forms an important wing of the army to provide firepower daring war 
as well as dining cross-border skirmishes with the enemy. Artillery generally falls into 
three basic categories; guns, howitzers and mortars. The principle difference between 
them being the trajectory of the round fired. A gun has a high muzzle velocity, and is 
normally used in a direct fire mode where the target can be seen and penetration is 
desirable. Howitzers have a somewhat lower muzzle velocity and arc their shells onto a 
target. They are used in both a direct fire and indirect fire mode. This is especially useful 
when an enemy is concealed behind a prepared position or the artillery men desire to 
have a shell explode over an enemy’s head. The airburst does less damage to hardened 
targets, but causes many more human casualties due to shell fingmentation covering a 
large area. Mortars have a very pronounced arc of flight. They have a relatively low 
muzzle velocity and are unsuitable for direct fire. Their principle value comes fi:om being 
able to lob shells behind an obstacle, such as a fortification or a hiU. They are not very 
accurate and dependent upon the amount of propelling powder to determine the point of 
impact. 

The effectiveness of artillery is largely judged by the accuracy in hitting the 
targets. The accuracy and reliability is influenced by the design criteria used in desigmng 
the shell. Artillery shells are a class of projectiles around which much of the aeroballistic 
theory was ori ginally developed, and they continue to form a significant part of the 
aeroballistician's interest. 
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Study of the motion of projectiles through an external medium is known by the 
name of external ballistics. If the external medium is the earth's atmosphere, then external 
ballistics become S5monymous with aeroballistics. 

The conventional approach hitherto for understanding the in-flight behaviour of 
projectiles was to develop mathematical models that could predict all elements of the 
trajectory from launch to target. For this purpose, it becomes essential that all forces, 
moments affecting the flight of the projectile be accounted for in a well-defined 
mathematical form^ Beginning with the most simple but relatively inaccurate 
mathematical model, the in-vacuo trajectory model, more and more sophisticated models 
of increasing accuracy such as the point mass model (PM), the modified point mass 
model (MPM) and the six-degrees-of-freedom (SDF) model have been developed^ These 
models require aerodynamic coefficients as input (e.g., drag coefficient , damping in 

roll derivative Ci^, lift curve slope etc.). Accurate values of these parameters are 

required in mathematical model for description of artillery shell dynamics. 

There are three distinct approaches for estimating aerodynamic parameters: 

• Theoretical methods 

• Wind tunnel testing 

• Flight-testing 

At the preliminary design stage of any system, theoretical methods are used in spite 
of their limited accuracy. Wind tunnel testing improves the accuracy of estimation but is 
a time consuming and expensive way of estimating aerodynamic parameters. Precise 
simulation of control surfaces, power effects of flight conditions is difficult. The model 
tested in the wind tunnel is generally slightly different from actual flight due to the last 
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minute configuration changes. Other reasons for discrepancies between flight and wind 
tunnel results are Reynolds number discrepancies and interferences due to support 
system. It is therefore desirable that the wind tunnel estimates be corroborated with the 
estimates from actual flight test data. 

The most commonly used methods for parameter estimation from flight data are 
broadly classified into the following categories. 

• Equation error methods 

• Output error methods 

• Filtering approach 

• Filter-error methods 

• Neural Network based methods. 

The principle of Least squares is used in the equation error method wherein the 
error gets minimized with respect to the unknown parameters in each of the equations. Its 
advantages include computational simplicity, non-iterative nature and applicability to 
both linear and non-linear models. The disadvantage is that the method caimot be directly 
applied if all the states are not measured accurately, and produces poor results if the 
measurements are noisy. Therefore in order to get accurate results, considerable effort is 
required for data reconstruction and smoothing. Sometimes these data processing tasks 
are more complicated than the task of parameter estimation itself. However these 
methods are useful as startup for more advanced estimators. 

In the output error method, the error between the measured and the model response 
is mi nimiz ed. The method assumes that there exist no modelling errors. This method 
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processes the measurement noise while assuming the model representation of the given 
system to be exact. A comprehensive survey of these methods is reported by Maine and 
Iliff* . The methods like analog matching^ , Newton Raphson method ’ , modified Newton 
R^hson method^ etc. fall in this category. Maximum likelihood (ML) method is one of 
the variants of output error methods. Maximum Likelihood estimates are those for which 
the observed value would be the most Ukely to occur; “most likely” is defined as 
maximization of likelihood function of the observed responses, given the parameters. The 
main advantage of this method is that parameter estimates are asymptotically unbiased, 
consistent and efficient provided the model assumed is correct. The schematic block 
diagram of ML method is presented below through Fig 1.1 and algorithm is presented in 
Appendix B. 



Fig. 1.1 Output error method 

Filtering approach*’’ is an extension of Kalman filter’®'" . Kalman filter is a linear, 
optimal estimator of state variables of a linear, time-varying system, operating in a 
Gaussian stochastic environment. Optimal estimator here is referred to a computational 
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algoritbm that processes measiirements to deduce a minim um error covariance of the 
state of a system combining all the information available i.e. 

° knowledge of system and measurement dynamics 
° assumed statistics of system noise and measurement 
° initial condition information. 

Traditional Kalman filter is applicable and optimal for linear systems only. When 
either the system or the messurement equations are non-linear, the same algorithm can 
still be applied by local linearization the system about the current state. Such filter 
applied to nonlinear systems is called Extended Kalman Filter (EKF)*’®’^^ , and it need not 
be optimal. For Kalman filter to be used as a parameter estimator, the unknown 
parameters are added as additional state variables, by considering the constant system 
parameters vector as the output of an auxiliary dynamic system. Parameter estimation is 
inevitably a nonlinear estimation problem whether the system and measurement 
equations are linear or non-linear. The detailed description about Extended Kalman Filter 
is presented in Appendix A. The schematic block diagram of EKF is presented below 
through Fig. 1.2. 


input 


process 

noise 


measurement 

noise 


Flight Vehicle 


■> 

■> 


Mathematical model 

Augmentation of states 
by unknown parameters 


Estimation of augmented 
states by extended 
kalman filter 
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measured 

response 


Fig. 1.2 Filtering approach 
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EKF produces estimates of the parameters that approximately mmimize the mean 
square error in the parameter estimates themselves as opposed to mi nimi zing a cost 
function that is based on matching the output variable behavior given a specific input 
trajectory, which is what most ML and least squares techniques do. This filtering 
approach is well suited for online applications, can handle process noise (modeling 
errors) and measurement noise (instrument errors). Extended Kalman Filter produces 
time histories for each of the estimated parameters, which is very useful if the time- 
varying nature of the parameters is of interest. But the performance of the filter strongly 
depends on the initial statistics of process noise, measurement noise and apriori 
covariance matrix. Information about the measurement noise covariance can be obtained 
from the laboratory calibration of measurement sensors used. Process noise is however 
more difficult to determine. Several procedures for obtaining measurement and process 
noise prior to estimation are available in literature. Well known methods by Mehra*^ and 
Morelli^'’ are being used with partial success. Other kind of approach known as adaptive 
filtering adjusts the covariances during estimation. Myers and Tapley’^ gave one such 
algorithm of adaptive filtering/tuning. Extended Kalman Filter is only a first order 
minimum variance filter. The higher order terms neglected in the propagation of the error 
covariance can possibly lead to biased estimates. Various techniques such as iterated 
EKF(IEKF) or second and higher order filters are often used in many applications to 
reduce the estimation errors. In the present work no attempt has been made to use either 
adaptive filtering or higher order filters. 

Filter error method*’ is the most general stochastic approach to parameter 
estimation, which accounts for both process and measurement noise. It can give good 
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estimates from data collected in turbulent atmosphere(process noise) also. In this direct 
approach an appropriately defined cost function, usually the statistically formulated 
likelihood flmction, is optimized with respect to the unknown parameters using some 
suitable optimization method. The Gauss Newton method is generally used for this 
purpose. Filter is used only for the natural purpose of obtaining the true state variables 
from the noisy measurements, i.e. for state estimation. The scheme describing this 
method is presented in Fig. 1.3. 



Fig. 1.3 Filter error method 


Application of Artificial Neural Networks (ANNs) for flight vehicle modeling and 
parameter estimation is a relatively new emerging field*. We have reported works on 
application to aircraft, not yet to rockets and artilleiy shells. In this approach the aircraft 
motion variables and control inputs are mapped to predict the total aerodynamic 
coefficients. Of the various architectures of NN, two types namely Feed Forward Neural 
network (FFNN) and Reccurent Neural network (RJNN) have found favour in parameter 
estimation. FFNNs lead to a black box type of modelling, where no physical significance 
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can be attached either to the network structure or to the weights. The weights or 
parameters of the network are estimated by the back propagation algorithm. On the other 
hand , RNNs are amenable to state space modeling. They have fixed number of mutually 
connected nodes equal to the number of unknown parameters of the postulated state 
space model. Partly due to the fixed structure and partly due to difficulties involved in 
tuning the sigmoidal nonlinearity of RNNs, FNNs found wider applications than RNNs. 

Parameter estimation of an aircraft with many measurements like, acceleration 
(both linear and angular), angular orientation, speed, angle of attack etc. has been 
discussed in detail in Ref.8. But firom cost effectiveness point of view (and also because 
of the non-availability of space for instrumentation), it may not be feasible to use many 
sensors for artillery shells going through many development flight trials. 

As an alternative approach, limited flight tests are routinely conducted with scale 
down model in Aero-ballistic range facility'^. The test facility is an enclosed, 
atmospheric, instrumented concrete structure used to examine the exterior ballistics of 
various free flight configurations. The facility contains gunroom, control room, model 
measurement room, blast chamber and instrumented range. The range may have large 
number of locations available as instrumentation sites. Each station is used to house fully 
instrumented orthogonal shadowgraph stations. Using this instrumentation, spatial 
positions, speed and angular orientations etc. are estimated. These are then used for 
parameter estimation. 

Winchebach et al'* conducted firee flight aeroballistic tests to obtain subsonic and 
supersonic aerodynamic parameters of a wrap around fin configuration. The aerodynamic 
parameters presented were extracted firom the position-attitude time histories of the 


8 



experimentally measured trajectories using non-linear numerical interpolation data 
reduction routines. 

Defense Research Establishment, Quebec, Canada is currently engaged in 
research on high length/diameter ratio energy penetrators for 70mm air launched rocket 
application. One of the objectives of the program is to enlarge the free-flight aerodynamic 
database of high l/d finned projectiles*’. Dupis investigated six model configurations 
(scale down) during the test program*’. Standard SDF routine with ML method was used 
to match the theoretical trajectory to the experimentally measured trajectory. In that 
analysis, as a first step, the basic dynamic range data (time, position, angles), physical 
properties and atmosphere conditions were assembled and then SDF (with ML) were 
performed to estimate the aerodynamic parameters. The results of this analysis showed 
that once time histories of position, angles are recorded it is possible to estimate 
aerodynamic parameters using free flight test data. 

It is worth mentioning that aerodynamic parameters are strong functions of Mach 
number. In order to capture this dependence on Mach number, projectiles are fired at 
different initial Mach number and the basic data corresponding to that Mach number are 
recorded and analyzed for parameter estimation. Though this approach looks attractive (if 
such facility is available), however it has the following two distinct limitations: 

1 ) It may not be possible to simulate actual flight path that the full-scale model would be 
traversing. 

2) It requires large number of firings to capture Mach number dependence of 
aerodynamic parameters. 
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In view of the above, there have been conscious efforts towards generation of 
flight data of free flight projectile by simultaneous measurements with pulse radar or 
optical system for position data and continuous wave Doppler radar for velocity 
measurements. In order to improve the overall accuracy several systems of each type 
have often been used in parallel. However, this complicated and ver\' expensive set up is 
not required today. From a single continuum wave Monopulse Azimuth and Elevation 
Doppler Tracking Radar it is possible to extract highly accurate velocity and position data 
in three dimensions without any assumptions of ballistic behavior. 

Defense Research Organization of India, is currently engaged in design and 
development of new artillery shell, rockets etc. There is a routine demand for post flight 
analysis of flight data. In particular, there is strong demand for estimation of aerodynamic 
parameters from flight data to validate the design as well as to prepare database for fire 
control system (FCS) to enhance artillery effectiveness. Accurate values of aerodynamic 
parameters are of paramount importance in building effective FCS. 

In the present work, ML and EKF method are used for parameter estimation 
purposes. The motivation here was to conduct a study on the applicability of these 
methods in extracting aerodynamic parameters by processing flight data obtained for 
different class of flight vehicle/store. These methods have been applied to flight data 
starting from a one-dimensional ballistic target to, an artillery shell, longitudinal motion 
of Short Range Strategic Missile (SM) and also to flight data of an aircraft. As discussed 
earlier, the flight data obtained through aircraft maneuvers contain large number of 
information regarding its motion and control variables. This is possible because 
exhaustive instrumentation for data acquisition is possible for such aircraft. Flight data of 
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one-dimensional ballistic target contain very limited information regarding its 
motion/control variables. Same is the case for flight data obtained through radar tracking 
the spatial coordinates and velocity of a shell/rocket in flight. Any parameter estimation 
method requires adequate information about vehicle dynamics to estimate aerodynamic 
parameters correctly. Unlike flight data of aircraft, others (rockets/shells) face this 
problem of lack of information content in relation to estimation of aerodynamic 
parameters. Thus challenges remain at place to evolve strategies to estimate aerodynamic 
parameters using such data of limited information content. This study systematically 
explores the possibility of applying the methods to estimate aerodynamic parameter 
under such constraints. Further, the flight data obtained through typical strategic missile 
vehicle has some flavour of both aircraft and rocket/shell. A strategic guided missile has 
control surfaces, lifting stufaces and is generally housed with more number of sensors 
like accelerometers, gyros etc. as compared to none for shell or rocket. Thus a typical 
flight data obtained through such flight vehicle have more information content as 
compared to free flight rockets or shells. However, excitation of such missiles with best 
control input (for the purpose of parameter estimation) may not always be feasible, where 
as it is always a routine task to generate flight data with such input for aircraft. Hence 
information content available with such missile data though superior in comparison to 
free flight rockets/shells but is many fold inferior as compared to aircraft flight data. 

The applicability of the existing methods (ML and EKF) is tested for these three 
different classes of flight data. It is observed that both ML and EKF can be 
advantageously applied to estimate few aerodynamic parameters with acceptable level of 
accuracy. 
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Chapter 2 describes methodologies used to generate flight data for type of vehicle 
considered. In chapter 3 results and discussions are presented and chapter 4 concludes 
with scope for future work. 
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CHAPTER 2 


GENERATION OF FLIGHT DATA 

Due to the non-availability of real flight data, simulated flight data are generated 
for the purpose of parameter estimation. Mathematical models used to generate flight 
data of ballistic target, aircraft longitudinal motion. Short Range Strategic Missile(SM) 
and artillery shell are presented in subsequent paragraphs. 

2.1 Ballistic target 

Knowledge of target ballistic coefficient can be used in advance guidance laws 
such as predictive guidance to relax the interceptor acceleration requirements. In 
addition, knowledge of the target ballistic coefficient is required for fire control due to 
the importance of accurate intercept point predictions in launching the interceptor on a 
collision course. Therefore, accurate estimation of ballistic coefficient of a target re- 
entering the atmosphere is very important for both guidance and fibre control purposes. 
The flight data for simulating such a vehicle motion is modelled to investigate the 
applicability of both ML and EKF method'® in extracting parameter (ballistic coeff) from 
flight data. 

To generate flight data, one dimensional motion has been considered*®. In this 
example it is assumed that a ballistic target is falling on a straight line path directly 
towards a surface-based tracking radar. Only drag and gravity force acting on the ballistic 
target have been modeled. 
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Figure 2.1 Forces acting on a one-dimensional ballistic target. 


In this problem, the tracking radar measures the distance from the radar to target every Ts 
seconds. Referring to Fig. 2.1, it can be seen that the drag acts upwards whereas gravity 
acts downward. The total acceleration acting on the ballistic target can be expressed in 
terms of a zero lift drag C^o or a ballistic coefficient, as 


^drag ref^ DoS 


( 2 . 1 ) 


where = 




, g is the acceleration of gravity and q is the dynamic pressure. The 


dynamic pressure can be expressed in terms of the air density p and ballistic target 
velocity 

^=0.5pFj?2 (2.2) 

For the purpose of simplicity, we can assxune that the air density varies exponentially 
with height and its variation is modelled using the following relation, 

(2.3) 

Therefore two differential equations that govern the one dimensional ballistic target can 
be summarized as 
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(2.4a) 


R =V 



^T2 

0.0034e ^^gVj .2 
2J3 


(2.4b) 


These equations were solved using 2nd order Runge-Kutta numerical integration 
technique. The relevant parameter (ballistic coefficient) used to generate the ballistic 
target flight data has been assumed to be 500 Ih! . Time histories of range ( Rj-j ), 
which will be referred as Ballistic Target Range data (FD-BT) for future reference and is 
presented in Fig. 2.2. This will be treated as measured flight data. This data will be 
processed using both ML and EKF method to estimate relevant aerodynamic parameter 
(ballistic coefficient). 



Figure 2.2 Time histories of Ballistic Range data corresponding to >9=5001b/ft^ . 
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2.2 Aircraft 


Knowledge of aerodynamic parameters are of paramount importance in 
postulating mathematical models for aircraft dynamics. Aircraft dynamics is generally 
dominated by two modes namely, longitudinal mode and lateral-directional mode. 
Longitudinal mode can have both short period(u=constant) and phugoid(a=constant) 
modes. In this exercise, only short period excitation is considered. For short period mode 
accurate values of at least eight aerodynamic parameters are required to describe its 
dynamics. The relevant aerodynamic parameters are C^,C^ ,Q ,0^^ 

and .In order to estimate these parameters, it is necessary to excite the airplane in its 

longitudinal mode by perturbing the aircraft with external input of a particular type. 
Elevator input could be of various types; pulse, double pulse or specifically designed 3-2- 
1-1 input. Excitation using 3-2-1- 1 input enhance information content in the flight 
data®’®. Typical elevator inputs used for simulating aircraft flight data (FD-A) are 
presented in Fig. 2.3(a) and Fig. 2.3(b). 


2j , , , , , ^ , 

1 ■ 

0 . 

.1 . 

*6 

-3 ■ 

-4 - 

-5 

■^0 2 4 6 8 10 12 14 16 18 20 

Ttme, sec 

Fig. 23(a) Elevator input typel 
Pulse input. 



Time, sec 


Fig. 2.3(b) Elevator input type2 
3-2-1-1 input. 
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The mathematical model used for simulating short period perturbed 


longitudinal mode is given below . 


pus 

2m 


2 /., " 


where the assumed aerodynamic model is given by 


(2.5) 

( 2 . 6 ) 


Q =6:^ +Q^.a + Ci^.(9c/2M) + Ci^.<J 

Short period equations of motion were solved using 4*'' order Runge-Kutta method. The 
geometric, inertia and flight conditions used to generate simulated data are given in Table 
2.1. The values of aerodynamic parameters used for generating simulated flight data are 
presented in Table 2.2. 

Table 2.1 Geometric/mass characteristics of aircraft 


SL.No. 

Geometric/mass characteristics 

1. 

Chord length 

: 5.6m. 

2. 

Wing surface area 

: 38.40 ni2. 


Mass 

; 9227.40 kg. 

3. 

ly 

: 67935.39 kg-m^. 

4. 

Speed 

: 61.72 m/s 


Table 2.2 Aerodynamic characteristics of aircraft. 



Q 

l-q 

c, 

Lg 

iriQ 



mg 

-0.0522 

2.6307 

4.4134 

1.0425 

-0.0163 

-1.16 

-1.1228 

-0.3557 
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Time histories of angle of attack and pitch rate corresponding to various elevator input 
are presented in Fig 2.4 and 2.5 and are termed as FD-ACl for elevator input type 1, and 
FD-AC2 for elevator input type 2. Time histories of a, ^ will be processed using ML and 
EKF method to estimate aerodynamic parameters. Further, the simulated flight data with 
varied degrees of measurement noise were also used to check the robusmess of the 
methods. Fig. 2.6 and 2.7 present graphically the FD-AC with 10% noise to be used for 
estimation purpose. 



Time, sec 


Figure 2.4 Simulated response a and q for step input, without noise (FD-ACl). 
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10 12 14 16 18 20 

le, sec 

‘or step input, with noise 10% (FD-ACl). 
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0.1 



Time, sec 


Figure 2.7 Simulated response a and q for 3-2-1-1 input, with noise 10% (FD-AC2). 


2.3 Short Range Strategic Missile (SM) 

Strategic Missiles are generally used for short range (2-4 km) engagements. It has 
three distinct phases of flight path. It has to climb to reach a particular height, cruise for 
certain distance and then finally hit the target under the commanded guidance and 
control. In this work the third phase is taken to be a ballistic trajectory i.e., without 
guidance. It is mitially propelled Avith high thrust, and subsequently after certain period of 
time, sustainer (with much lower thrust) takes over for a desired period of time. However, 
during the initiation of guidance and control, the missile is without any active thrust. In 
order to achieve predefined flight path, a set of preset control inputs are applied. The 
dynamics of this motion has been modeled using SDF trajectory model. In this study only 
the longitudinal motion is being considered. The equations representing SDF model are 
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listed in the subsequent section. The mass, moment and aerodynamic characteristics used 
to solve these equations of motion are listed in Table 2.3 and Table 2.4 respectively. 


Table 2.3 Geometric/mass characteristics of Strategic Missile. 


SL.No. 

Geometric/mass characteristics 


Length/Dia 

: 8.33. 


Mass, initial/final 

: 1.258. 


lx, initial/final 

: 1.2325 


ly, initial/final 

: 1.1657. 


Iz, initial/final 

: 1.1657. 


Table 2.4 Aerodynamic characteristics of Strategic Missile. 




Cr 

nig 

Cm 

mg 

0.3 

5.0 

2.35 

-6.5 

-8.9 


The Six degrees of freedom model used to represent the dynamics of SM is 
modeled using the following sets of equations as presented below. 


u= {^s! m)C^ — qw ■¥ rv - gsmO + {^l ni) 

(2.7) 

v= {iqsl Tn)Cy - ru + pw + gsin<p cos6 

(2.8) 

w= (qs/m)C^ - pv + qu + gcos<j^ cosO 

(2.9) 

1 

+ 

fl 

(2.10) 

q={qsdC„+rp{l^-I,)}/ly 

(2.11) 

r=\^qsdC„+pq{l,-Iy)}/l, 

(2.12) 
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/> + g tan0 sin^+ r tan^ cos<^ 
qcos(f> — rsm.^ 

\j/= rcos^zisec^ + ^sin^isec^ 


(2.13) 

(2.14) 

(2.15) 


To derive the spatial position equations, above equations were used to transform 
the body-axis velocity (m, v, w) into earth fixed axis. The equations are: 


Jr= wcos^4^cos^ + V (cos^^/■sin^sin^-sin^^^cos^) 
+ w sin^cos^+sin^ sin^) 

T=MSuni/'cos^ + V (sin^^sin^sin<^ + cos^i/^cos(Z>) 
+ w {siny/sinOcos^-cosyrsin^) +W^ 

Z=usin0 — V cosffsin^ — wcosffcos^ +W, 


(2.16) 

(2.17) 

(2.18) 


Aerodynamic model used in this analysis is as given below: 

cosacosyJ+Q sin a cos a sin/? 

^y=Cyy.ind smp 

Cz- -Qcosa -C^^sinacos/?-C^^sinasin;6 

where, 

Q = Qo +q„a+C,, {qdllV) + 

^ywind C„+C,fP*C„(rdl2V) 

^Dwind ^d 

C,=c„ (pd/iv) 

’^^ma l2V^ + 

C,=C^+C,,P^-C„(rdl2r) 


Wind model used in this analysis is as given below: 

w'= u-W^ COS0 cosy/- cosO smy/ sin^ 
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v'—v—W^ (cos^i^sin^sin^ii— siiny cos^) 
-Wy (cos^^/^cos^+siin«^sm^sin^) 
-W^cosOsin^ 

w'=w—W^ (cos^sin^cos(^ + siin4Asin(!>) 
-Wy (siini/-sm0cos^-cos^^/^sin^) 
-W^cos9cos(f) 


thus. 


a=tan"’ 



and 


;5=tan ‘ 



where are the wind components blowing toward the 


X, y, and z directions. 


The flight data so generated corresponding to the elevator input as given in Table 2.5 are 
pictorially presented in Fig. 2.8. This flight data will be referred as FD-SM while 
estimating parameters from flight data representing flight of a short range strategic 
missile. 



Time, sec 


Figure 2.S Measured variables with respect to time without noise (FD-SM). 
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Table 2.5 Control surface deflection with time. 


Time interval 

Control surface 

deflection, degrees 

0<t<8.5 

-10.6 

8.5 <t< 16.25 

10.0 

12 <t< 16.25 

-4.5 

16.25 <t 

0.0 


2.4 Artillery shell 

In general, the trajectory followed by center of mass of artillery shell in space is 
three-Dimensional {x,y,z). Theoretical modelling of such trajectory demands 
aerodynamic modeling to be exhaustive. It requires usage of large number of 
aerodynamic parameters. Hence estimation of all parameters at one time might be 
involved. Initially in estimation algorithm only Point Mass Model (PM) was considered 
and as the work progressed , the trajectory model was made 3 -Dimensional using 
Modified Point Mass Model (MPM) and Six Degree of Freedom model (SDF) 
respectively. 

2.4.1 Point Mass Model (PM) 

In this model, it is assumed that the only aerodynamic force acting on the 
projectile is drag. It provides fairly accurate estimates of range for adequately stable 
projectiles and can also be used to estimate the first order effects of wind. The equations 
of motion for point mass model are given as below: 


24 




fx_ Tvpd^C^ ^/ 
dt ^ 8/m 

]2 


dt 


d y 




dP '' g/M 


dP 


8/m 


dt 


-w. 


(2.19) 


I 


^ dy 
[~di~ 

\ 

w 

y 5 

J 

(2.20) 

\ 




• 

(2.21) 


where x denotes range, y denotes height, z denotes drift and W^,Wy, respectively, 
the X , y and z components of the wind velocity W . This model, however, does not 
account for the effect of spin of the shell and thus fails to predict the drift due to yaw of 
repose moment. Also, it neglects the lift forces acting on the shell. 


2.4.2 Modified Point Mass Model (MPM) 

The modified point mass model is also known as four degree-of-fireedom model 
(three spatial degrees-of-fi:eedom plus axial spin). Its basis is a conventional point mass 
model, in addition, the instantaneous equilibrium yaw is calculated at each time step 
along the trajectory so as to provide estimates of yaw, drag, and drift resulting from the 
yaw of repose. The equations of motion for modified point mass model are given as 
below* : 


du 

~dt 


Ttp d^ 
8/m 




+ 


Ttp d^ 

8/m 




- g’o 




dp _ Ttp d* 

~dt~ 161, 




(Drag) 

(Lift) 

(Gravity and Coriolis) (2.22) 
(Spin damping) (2.23) 
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(Yaw of repose) 


(2.24) 


-Zpl^ \vx[duldt) 

ff = 1 i i= i iri 

^ npd^C^, 


The quantities in the above equations are defined as follows: 

— dx ~ ~ - - -r 

u = — , V=u—W, r=x—R 
dt 

if =(0,-i?,0), [IF Denotes the wind vector] 

CO =(-Q cos[latitude] cos[azunuth], -Q sin [latitude], Q cos[latitude] sin[aziniuth] ) 

where, Q = 7.29x10'^ rad/s. (rotation of the earth), 

R = 6370320 m (radius of the earth), 

go = 9.80665 [1 - 0.0026373 cos (2xlatitude) + 0.0000059 [cos (2x latitude)]^], 
the axis system is as for the point mass model with x being along the line of fire. 

As may be seen firom the above equations, the aerodynamic coefficient input 
required is extensive and accuracy of these aerodynamic coefficients is crucial for 
reliable estimates of range. However, for most artillery shells, the reliability of estimated 
values of aerodynamic coefficients is not high enough to inspire confidence in resulting 
range estimates firom this model. In spite of this lumtation, this model has the capability 
of providing reasonable accounts of end point data. It can estimate wind corrections and 
can accept non-linear aerodynaimc inputs if required. 


2.4.3 Six Degree of Freedom model (SDF) 

'fhe most sophisticated trajectory model developed was the sk-degree-of-fi:eedom 
model having the spatial degrees of freedom, yaw degree of fi'eedom in two planes and 
the spin. In this model there are no assumptions concerning linearized aerodynamics or 
projectile symmetry. However, the indeterminability of many of the uutial conditions and 
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aerodynamic coefficients which are required as input frequently results in the model not 
giving sigmficantly better end results than the modified point mass model. Thus its usage 
might not be justified for routine fire control work. It is nevertheless a powerful tool for 
the ammumtion designer. The equations of motion for six-degree-of-freedom model are 
as given for Strategic Missile except a small change in aerodynamic modeling. Here the 
aerodynamic modeling doesn’t include the aerodynamic derivatives due to control 
surface deflection S as there are no control surfaces attached to shell. The description of 
SDF model has already been presented in section 2.3. The mass and aerodynamic 
characteristics of the artillery shell used in generating flight data are Table 2.6 and 2.7 
respectively. 


Table 2.6 Geometric characteristics of the artillery shell. 


SL. No. 

Geometric characteristics 

1. 

Total length 

: 825 mm. 

2. 

CG location 

: 533 mm from nose tip. 

3. 

Diameter 

: 155 mm. 

4. 

Mass 

: 42.6 kg. 

5. 

Moment of Inertia 


Axial, 4c 

: 0.146 kgW . 


Traverse, 4y 

:1.709kg-m^ 
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Table 2.7 Aerodynamic characteristics of the artillery shell. 


SL. No. 

Mach 

Number 


Qa 



1. 


0.13 

3.2 

mHi 


2. 

0.85 

0.13 

3.2 

ngunm 


3. 

0.9 

0.13 

3.25 

1.2 

-0.03 

4. 

0.92 

0.14 

3.3 

1.1 , 

-0.03 

5. 

0.96 

0.203 

3.45 

0.95 

-0.03 

6. 

0.98 

0.269 

3.5 

0.95 

-0.03 

7. 

1 

0.329 

3.5 

1 


8. 


0.351 

3.4 

1.1 

-0.0295 

9. 

1.05 

0.367 

3.25 

1.55 

-0.029 

10. 

1.1 

0.366 

3.15 

1.85 

-0.029 

11. 

1.15 

0.353 

3.15 

1.95 

-0.0285 

12. 

1.2 

0.343 

3.14 

2 

-0.0282 

13. 

1.3 

0.326 

3.13 

2.07 

-0.028 

14. 

1.4 

0.314 

3.12 

2.15 

-0.027 

15. 

1.6 

0.293 

2J5 

2.2 

-0.026 

16. 

1.8 

0.273 

2.5 

2.3 

-0.025 

17. 

2 

0.257 

2.48 

2.32 

-0.024 

18. 

2.2 

0.245 

2A7 

2.35 

-0.023 

19. 

2.4 

0.236 

2.45 

2.35 

-0.022 

20. 

2.6 

0.229 

2.4 

2.35 

-0.0215 


It may be mentioned that the flight data used for parameter estimation (for shell) 
was simulated using SDF model. A typical flight data showing X,Y,Z,u,v,w,p,q,and r has 
been presented in Fig. 2.9. However during estimation depending upon the number of 
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aerodynamic parameters to be estimated, selectively PM, MPM, SDF models were used. 
The flight data for artillery shells will be referred by FD-SDF for future reference. 



Fig. 2.9a Typical measured data of X, Y, Z without noise. 





Fig. 2.9b Typical measured data of u, v, w without noise. 
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CHAPTERS 


RESULTS AND DISCUSSIONS 

In this chapter, estimated parameters via ML as well as EKF methods are 
presented. The simulated flight data corresponding to Ballistic target (FD-BT), Aircraft 
Strategic Missile (FD-SM) and Artillery Shell (FD-SDF).are used measured 
flight data. The effect of measurement noise and initial variance on estimated parameters 
is also studied. The estimated parameters are presented only with its standard deviation to 
asses the accuracy of the estimates. 

3.1 Estimation of ballistic coefficient from ballistic target data (FD-BT) 

As stated earlier, here the ballistic target was assumed to fall on a straight-line path 

directly towards a surface based tracking radar. The ballistic coefficient ( B - — — — ) is 

^ref^DO 

estimated using both ML and EKF method. Table 3.1 compares the estimates obtained 
using both the methods. 

Table 3.1 Comparison of estimates of ballistic coefficient obtained using both ML 
and EKF. 


Initial guess 
value 

Measurement 

Noise 

variance 

Estimated 
(BaUistic coel 

parameter 
fiicient) using 

ML 

EKF 

600 

500 

500.1330 

(0.006163)* 

498.2207 

(0.070775)” 

600 

1000 

500.1881 

(0.006164) 

498.1900 

(0.070576) 

1000 

500 

500.1330 

(0.006163) 

497.7265 

(0.069323) 

1000 

1000 

500.1881 

(0.006164) 

497.3556 

(0.068497) 


*CR-Rflo bound. 

° Standard deviation. 
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During estimation, the initial guess values were varied from 600 to 1000 (lb/ft2) 
and noise variance were also introduced (500 to 1000 ft^ ). It can be well observed that 
both the methods estimate the value of Ballistic coefScient very close to the actual value 
of Ballistic coefficient (500 lb! ft^) used for generating the simulated data. Further both 
the estimates have low standard deviation showing higher confidence of the estimates. It 
is also noted that, for the particular case estimates are not dependent on the value of 
initial guess values. Further, both the methods show robustness with respect to noise 
variance. 

3.2 Aircraft flight data (FD-AC) 

Knowledge of aerodynamic parameters is of paramount importance, to develop 
accurate mathematical model representing longitudinal dynamics (short-period). The 
aerodynamic parameters needed to develop the aerodynamic model are force derivatives 
and moment derivatives and . Thus, a study was 

carried out to explore the possibility of extracting the parameters using both ML and EKF 
method. Researchers have already presented applicability of both these methods to such 
flight data*-’ .However, in this study, this case is included to validate the in-house 
developed code. Further the effect of noise and type of input on the estimates were also 

studied to check the robustness of the method . 

The flight data (a, q) FD-ACl and FD-AC2 (as given in Fig. 2.6 and Fig. 2.7) 
were used to estimate the aerodynamic parameters by both ML and EKF method for 
varied degrees of variance and percentage noise (5% and 10% being presented). Table 
3.2 presents the numerical values of the estimated parameters with standard deviations. 
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Table 3.2 Comparison of aerodynamic parameters of obtained by using ML and 
EKF method on FD-ACl and FD-AC2 with 5% 
measurement noise. 


Parameters 

True 
values of 
parameters 

FD-ACl 

FD-AC2 

EKF 

MMLE 

EKF 

MMLE 


2.6307 

2.6324 

2.6681 

2.6170 

2.6239 



(0.146776) 

(0.1317) 

(0.090307) 

(0.04294) 


4.4134 

4.2655 

4.3077 

4.3591 

4.2334 



(1.005723) 

(1.214945) 

(0.672155) 

(0.367366) 


1.0425 

1.0436 

1.0556 

1.0367 

1.0396 



(0.04696) 

(0.075803) 

(0.028275) 

(0.036135) 

C 

-1.16 

-1.1588 

-1.1593 

-1.1588 

-1.1581 

"‘a 


(0.008557) 

(0.009838) 

(0.00569) 

(0.002952) 


-1.1228 

-1.1122 

-1.1019 

-1.1254 

-1.1267 



(0.069582) 

(0.079352) 

(0.042872) 

(0.024079) 


-0.3557 

-0.3541 

(0.002741) 

-0.3549 

(0.007003) 

-0.3548 

(0.001829) 

-0.3551 

(0.002652) 


-0.0522 

-0.0520 

(0.002203) 

-0.0515 

(0.003092) 

-0.0525 

(0.002039) 

-0.0525 

(0.003094) 

^ Wo 

-0.0163 

-0.0162 

(0.000134) 

-0.0162 

(0.000478) 

-0.0163 

(0.000125) 

-0.0162 

(0.000479) 


Referring Table 3.2, it can be observed that both ML and EKF successfully 
estimated all the parameters. It can be seen that consistently, that estimated parameters 
(both ML and EKF') using 3-2-1-1 input yields lower values of standard deviation 
suggesting more confidence in these estimates. Similar observations could be extended 
for estimated parameters obtained processing flight data with higher degree of noise 
( 1 0%) as presented in Table 3.3. 
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Table 3.3 Comparison of aerodynamic parameters obtained by using ML and EKF 
method on FD-ACl and FD-AC2 with 10% 
measurement noise. 


Parameters 

True 
values of 
parameters 

FD-ACl 

FD-AC2 

EKF 

MMLE 

EKF 

MMLE 


2.6307 

2.5248 

(0.284995) 

2.7294 

(0.132548) 

2.6771 

(0.179135) 

2.6546 

(0.04305) 


4.4134 

4.6236 

(1.787376) 

4.4987 

(1.227096) 

3.8604 

(1.274208) 

3.4836 

(0.370112) 


1.0425 

1.0331 

(0.091799) 

1.0774 

(0.075854) 

1.0644 

(0.056584) 

1.0462 

(0.036349) 


-1.16 

-1.1613 

(0.015394) 

-1.1584 

(0.009927) 

-1.1548 

(0.010797) 

-1.1529 

(0.002944) 


-1.1228 

-1.1719 

(0.135488) 

-1.0874 

(0.080165) 

-1.1049 

(0.085014) 

-1.1138 

(0.024137) 


-0.3557 

-0.3567 

(0.005102) 

-0.3545 

(0.007017) 

-0.3553 

(0.003588) 

-0.3552 

(0.002653) 

Cl 

M) 

-0.0522 

-0.0525 

(0.004255) 

-0.0501 

(0.003089) 

-0.0527 

(0.004049) 

-0.0534 

(0.003111) 


-0.0163 

-0.0163 

(0.00025) 

-0.0162 

(0.000477) 

-0.0163 

(0.000245) 

-0.0162 

(0.000477) 


The convergence of the estimated parameters and variance with EKF is pictorially 
presented in Fig. 3. Land 3.2. It can be seen that after around 2-5 sec, the error variance 
converges to zero suggesting accurate estimates. 
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3.3 Estiination using flight duts of a typical short range strategic missile 

The estimation of aerodynamic parameters from flight data of such a missile is 
more involved as compared to that of an aircraft. It is not always possible to excite the 
missile with a best input (3-2-1-1). Further due to limitation of having large number of 
sensors for data acquisition (dedicated for parameter estimation), there is an inadvertent 
compromise on the quality and quantity of information content in such flight data. It is 
therefore, important to identify the minimum number of parameters required to be 
measured to initiate parameter estimation. Due to lack of information content, many a 
times the estimated parameters tend to become sensitive to the initial guess values. Thus 
a systematic study was carried out to identify the most suitable combinations of minimum 
numbers of motion/control variables required for parameter estimation. Various 
combinations of X, Y, Z, q,(t, and a. were selected and parameters were estimated 

using various values of initial guess values as presented in Table 3.4. Based on this table, 
it was decided to have X, Z and a. as motion variables to be used for estimation using 
ML method. Such exercise was also carried out for EKF method and decided to have X, 
Z and as the most appropriate variables for this purpose. It was further observed that the 
effect of pitch rate {q) on X and Z were too insignificant. Thus the estimation of from 

these sets of data was not carried out. Such a study using simulated data would be useful 
in designing the experiment for acquiring the relevant parameters for parameter 
estimation. However, it may be stressed that such observations must be validated using 
real flight data. For estimation purpose, ML method used flight data (X, Z,a, ) as 
presented in Fig. 2.8. Estimated values of aerodynamic parameters with standard 
deviation are presented in Table 3.5. For both cases of noise level (5% and 10%), ML 
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method yields lower values of standard deviation as compared to EKF. A study was also 
carried out to see the effect of initial guess values on the estimated aerodynamic 
parameters. Typical sets of initial guess values (as presented in Table 3.5) for 
aerodynamic parameters were chosen. It was found that ML method failed to estimate 
parameters for the initial guess values-2. However EKF could estimate even with such a 
deviation in initial guess values. Thus it is observed that for data with limited 
information, EKF may be used advantageously. 

Table 3.4 Range of initial values for which ML method was able to estimate 
aerodynamic parameters of a Short range strategic missile. 


Combination 
of variables 

Range of initial values 
expressed in percentage. 

a; 2 


X,Z.q 

-15,12 


-15,18 

X, Z,a. 

-31,37 

X,Z,a 

-12,13 

ii, .Cl, .a 

-10,13 


-12,10 

a, rO, q 

-13,13 

a. ,a, q 

-10,10 

a^,a,,a,q 

-10,13 

X, Z, a, , a. 

-21,38 

Tz a,, a,, a 

-21,38 

X, Z. a^,a^,q 

-21,38 

\ X. /. q. a 

-11,19 
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Table 3.5 Comparison of estimated parameters (SM) obtained using ML and EKF 
method for different initial guess values and measurement noise of 5%. 


Parameters 

True 

value 

Initial guess values-1 

Initial guess values-2 

MMLE 

EKF 

MMLE 

EKF 


0.3 

0.3000 

0.2999 

Could not 

0.29984 



(5.00E-10) 

(1.85E-05) 

estimate 

(4.93E-05) 


5 

4.9968 

5.0321 


5.0587 



(4.81E-05) 

(0.024639) 


(0.068510) 

Cl 

2.35 

2.3459 

2.3967 


2.4839 



(1.05E-04) 

(0.034804) 


(0.096597) 

c„ 

-6.5 

-6.4998 

-6.5012 


-6.5042 



(9.30E-08) 

(0.001081) 


(0.003042) 

c^ 

-8.9 

-8.9000 

-8.9064 


-8.9660 

”*0 


(3.56E-06) 

(0.006757) 


(0.010062) 


Parameters C, C, 

Initial guess values of parameters 1 - 0.2, 3, 3, 
Initial guess values of parameters 2 - 0.2, 2, 4, 


-5, -7 
-3, -4 
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3.4 Estimation of aerodynamic parameters using flight data of a typical 
artillery shell 

Due to the non-availability of radar tracked data, simulated flight Ha t ? ) were 
generated using SDF model. Trajectories were generated corresponding to dilferent 
launch angles. During estimation, the algorithm included PMM, MPMM and SDF 
progressively. Point mass model can be used advantageously to estimate drag coefficient 
(Cj); however to estimate other derivatives ), the estimation algorithm 

should include higher order of trajectory models, namely MPM or SDF model The 
effects of process and measurement noise in the estimation of aerodynamic parameters 
using EKF method has been presented in detail and are explained case by case. 

3.4.1 Estimation using Point Mass Model: 

Case 1: 

During simulation, the drag coefficient along with other parameters was kept as 
fixed numbers. Their variation with Mach number was not modeled in simulating the 
trajectory using SDF model. Constant values of aerodynamic parameters used are given 
in Table 3.6. 


Table 3.6 Values of aerodynamic parameters used for simulation in case 1 






0.3 

2.0 

3.0 

-0.025 


It may be noted that constant approximation for aerodynamic parameters are justified for 
subsonic flight regime. The estimation algorithm used only Point mass Model and hence 
expected to estimate the drag coefficient, Q . Range (X) and height (Y) corresponding to 
launch velocity of 818 m/sec and elevation of 7.6 degrees were used as measured 
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variables. During the implementation of EKF method, numerical differentiation was not 
used to linearize the state equations; instead the linearized matrix was computed using 
differentiated terms obtained analytically. For the computation of state transition matrix, 
first three terms were found to be sufficient. Initially the algorithm was run with 
measurement noise corresponding to 20m randomly distributed over the whole range 
without any process noise (assuming the model to be perfect). The estimation process 
converged to yield a value of 0.28858 for and convergence of is graphically 
presented in Fig. 3.3. 



Fig. 33 Convergence of Cd as plotted w.r.t. range , easel. 
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Fig. 3.4a Comparison of measured and 
EKF predicted trajectories easel. 



Fig. 3.4b Error in range and height 
for easel. 


The estimated trajectory obtained during EKF run is compared with the measured 
trajectory as shown in Fig. 3.4a. Further, the growth of error in range and height is 
presented in Fig. 3.4b. It can be observed that, the maximum error in range is about 50m 
where as in the height about 10m. The estimated value of Q is 0.28858 instead of 0.3. In 

order to improve trajectory matching and the value of estimated Cj, pseudo-process 

noise with variance 0.05 and 0.5 were added to X and Z respectively and estimation 
algorithm was re-run. It was interesting to observe that the matching between measured 
and estimated trajectories improved considerably, however the estimated Cj value 

deteriorated significantly. The convergence of , error m range and error in height as 

presented with respect to range presented in Fig. 3.5, 3.6 respectively. It can be observed 
that the maximum error has come down to 10m in range and 2m in height. However, it 
can be seen from Fig. 3.5, that the value converged to 0.2369, instead of actual value 
of 0.3 showing further deterioration. 
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Fig. 3.6a Comparison of measured and Fig. 3.6b Error in range and height 

EKF predicted trajectories. 


In order to improve value, a pseudo-noise of 1 x 10"^ was also added to and 


the flight data X and Z were processed to estimate the drag coefficient Cj . In this run, 
process noise for X and Z were kept fixed at 0.01 and 0.1 respectively. The estimated 
value of showed remarkable improvement by converging to actual value of Cj used 

for simulation. This improved convergence is presented pictorially in Fig. 3.7. Further it 
was observed that the maximum error in range and height remained well with in 10m and 
2m respectively. 
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Fig. 3.7 Convergence of Cd 

Case 2: Estimation of Cj V/s Mach number using Point Mass Model: 

The artillery shells when fired are expected to travel Mach number covering subsonic to 
supersonic regime. The aerodynamic parameters are strong functions of mach numbers. 
Therefore, it is important to extract aerodynamic parameters as a function of Mach 
number. In this case EKF was applied using the PMM to estimate functional dependence 
of Cj . The simulated data were generated by using aerodynamic parameters (as function 

of Mach number) given in Table 2.7. Process noise was added to X, Z and (0.01, 0.1, 
0.0001) respectively. EKF routine was run to obtain estimated value of range, height and 
Cj as function of Mach number as presented in Fig. 3.8 and 3.9 respectively. 
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Fig. 3.8b Error in height. 



Referring Fig. 3.8(a) and (b), it can be seen that the estimated max error in range and 
height are well with in 7m and 1.5m respectively. The estimated value of as function 
of Mach number presented in Fig. 3.9 compares well with the manufacturer supplied 
value of Cj Vs Mach number. Thus, it can be concluded that EKF can be advantageously 
applied to extract functional dependence of drag coefficient on Mach number. Similar 





observations, related to application of ML method to the above mentioned cases have 
been reported in ref (20). 

3.4.2 Estimation of aerodynamic parameters using MPM: 

Due to non-availability of measured flight data, simulated flight data (FD-SDF) 
were first constructed using constant aerodynamic parameters and then using 
aerodynamic parameters as function of Mach niunber. However, in the estimation 
algorithm MPM model was used for this case. 

3.4.2. 1 Flight data simulated with constant aerodynamic parameters: 

Aim of this study was to investigate the capability of EBCF method in extracting 
constant aerodynamic parameters using measured trajectory coordinates X, Y and Z. 
Successful implementation of ML method ( in extracting the above 

parameters are reported in Ref 20. However during estimation has been kept fixed to 
its true value. The EKF method was implemented and was successful only in estimating 
the parameters correctly. However the parameters and showed poor 

convergence. In order to improve their estimates, initial value for was set at the exact 
value, and it was seen as presented in Fig. 3.10. that the convergence of Cj^ improves 
remarkably. 

3.4.2.2 Flight data simulated using aerodynamic parameters varying with 
Mach number. 

The concept of using aerodynamic parameters as constant values is strictly valid 
for low subsonic flight. However the shells under investigation are supposed to cover 
flight regime covering subsonic to supersonic flow. Thus it is important to capture the 
Mach number dependence of these parameters in developing accurate mathematical 
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model. During implementation of EKF method, it was observed that it was not possible to 
extract all these parameters simultaneously using only X, Y and Z. Several trials were 
attempted to arrive at a particular combination of motion variables, which could be used 



Range, m 


Fig. 3.10 Estimation of constant parameters using MPM from FD-SDF. 
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Fig. 3.11 Estimation of C^and C,p as a function of Mach number 

(45 degrees elevation). 

to estimate only few aerodynamic parameters as function of Mach number. The flight 
data used for this purpose included X, Y, Z, u, v, w and p. The implementation of EKF 
resulted in estimation of and as function of Mach number, the comparison of 

C^and Cip estimated through EKF and actual are presented in Fig 3.11. However the 
mach number dependence of and could not be captured with this flight data. One 
possible reason could be due to the fact that the trajectoiy angle of attack of shells is well 
with in 0-0.3 degree for most of the flight path. Thus, the effect of and on range 

and deviations are too small to be sensitive for their estimation. There is a need to 
carryout detail study to modify the strategy to capture these variations. For the sake of 
completion, it may be mentioned that similar observations were noted while 
implementing EKF using SDF model. 
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During the estimation process, it was observed that ML methods are hi ghl y 
sensitive to the correctness of initial values of motion variables, however, EKF can 
advantageously applied by appropriately modifying process noise to handle such inherent 
error present in acquiring trajectory data. This is perhaps a imique advantage, with EKF 
which designers utilized rigorously for implementation of on-line control law and 
guidance schemes. 
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CHAPTER 4 


CONCLUSIONS AND SUGGESTIONS FOR FUTURE WORK 

4.1 Conclusions 

In the present work, ML and EKF methods have been applied to estimate 
aerodynamic parameters from simulated radar tracked data. These methods have been 
applied starting from a one-dimensional Ballistic target to an artillery shell, longitudinal 
motion of strategic missile and also to flight data of an aircraft. It is observed that both 
EKF and ML methods can be applied successfully to estimate all the parameters from 
flight data of an aircraft. However for shells/rockets, both the methods are capable of 
estimating only drag coefficient, and roll damping coefficient accurately. Both the 

methods failed to estimate the aerodynamic parameters and . However, for SM, 

it was found that, if appropriate flight data is generated, both the methods can be used 
advantageously to estimate force and moment derivatives. 

4.2 Suggestions for future work 

1. In the present work, simulated flight data were used to estimate aerodynamic 
parameters of the artilleiy shell and strategic missile due to non-availability of the real 
flight data. The capability of these methods should however be validated using real flight 
data. 

2. Scope of EKF may be expanded to capture Mach number dependence of all the 
parameters. 
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3. The simulated data was generated for standard atmospheric conditions. Applicability 
of these methods to estimate parameters using flight data under non-standard atmospheric 
conditions may be explored. 

4. Filter error method which uses Kalman Filter for state estimation and ML method for 
parameter estimation could also be explored. 


A...L45L4S1 — 
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APPENDIX A 


EXTENDED KALMAN FILTER 

Kalman filter applied to non-linear systems is known as EKF. Parameter 
estimation is inherently a non-linear problem whether the system is linear or non-linear. 

The dynamic system, whose parameters are to be estimated, is assumed to be 
described by the following stochastic equations. 

^(0 = /[^( O , w ( 0 ] + HO ^(^0 ) = ^0 

z(A) = /7[x(k),u(k)] + v(k) k=l,....N 

First equation represents the system dynamics and the second the measurement 
dynamics. The n and m dimensional system functions/ and h are general nonlinear real 
valued vector functions. These system functions are assumed to have sufficient 
differentiability to be able to invoke Taylor series expansion. The m-dimensional 
measurement z is sampled at n discrete time points with a uniform sampling time At. 

It is assumed that the process and measurement noise w and v affect the dynamic system 
linearly. They are assumed to be characterized by a zero-mean white Gaussian noise with 
covariances of Q and R respectively. 

The above dynamic model is for state estimation. For parameter estimation the 
unknown parameters 0 are added as additional state variables considering them as the 
output of an auxiliary dynamic system 

^ = 0 

The augmented state vector X is defined as = {x^ ,0^ } . 
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Now the extended system can be represented as 




f[x{t),u(t)\ 

1 - 

w{t) 

0 


_0 


z(k) = h^[X(k),u{k)] + v(k) 

where the subscript a denotes augmented system vectors. 

EKF algorithm: 

The initial conditions X(0), P(0) and the noise covariance matrices Q and R need 
to be known before applying the algorithm. In the present theses process noise is assumed 
to be zero and the simulated data is generated with constant measurement noise R. 

Step 1: Time update (effect of dynamics) 

= Xk* + ]f[X(t),u(t)¥t 

h-\ 

~ ^(k,k-\) ^k ^{k,k-\) + Q 
Step 2: Measurement update (effect of measurement) 

K, =P-HI{H,P;HI +R] 


Superscripts and ‘+’ denote the predicted and corrected state 

estimates/covariances respectively. Step 1 and step2 are also called as prediction and 
correction steps in Kalman filter. 
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The state transition matrix is given by e^' 

where F is given by 


also ^ k 


H, = 


dux{k\u{k)] 

dX{k) 

dK{X{k\u(k)\ 

dX{k) 


x,=xu 


X~ is the predicted state vector from integration over the time step of measurement 
data 

X* corrected state vector using the measurements. 

P~ predicted/apriori error covariance matrix. It gives error in the predicted states. 

P* corrected/posteriori error covariance matrix. It gives the error in the corrected 
estimates. 

R the measurement noise covariance, determines how much information from the 
sample is used 

Q the model/input noise covariance contribute to the overall uncertainty of the estimate 
as it is added to P~ in each time step. 

K the Kalman gain determines how much of the innovation (the difference between the 
actual measurement and model measurements) is used to correct the estimate. 
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APPENDIX B 


MAXIMUM LIKELIHOOD (ML) METHOD 

The parameter estimation algorithm used in this study is the conventional 
‘Maximum Likelihood Method’. The programs were written in ‘C’ using the MLE 
algorithm. The algorithm is given as follows: 

Let a mathematical model of a system be represented by the following equations: 

X(t,) = AX(t,) + BU(l,) (A.1) 

Z(t,) = CX(t,) + 0>,, (A.2) 

where X is the state vector, U the control vector, Z the measured (observation) 

response and is the noise vector respectively. The matrices A and 5 contain the 

unknown stability and control derivatives (parameters). If there is no state noise and the 

matrix GG is known, then the Maximum Likelihood estimator minimizes the cost 

function, 

y(0) =(1/2) {[ z((,)-z,(r,) J (GG^J' [ ]} (A.3) 

/=1 ^ 


where GG^ is the measurement noise covariance matrix and Zg{t, ) is the computed 
response estimate of Z at for a given value of the unknown parameter vector 6 . The 
cost function is a function of the difference between the measured and computed time 
histories. In our case, we have taken GG^ = I for the cost function. To minimize the cost 
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function J{6), we can apply the Newton-Raphson algorithm, which chooses successive 
estimates of the vector 0 of unknown coefficients. If k is the iteration number, then 


Jfe ) ]"' [ vj J(e, ) ] , (A.4) 

where the first gradient is defined by 

(A.5) 

and the second gradient is approximated by 

(GG’')''[V,Z,(r,)]}. (A.6) 


The matrix GG^ is approximated by a diagonal matrix, with the diagonal 
elements given by, 


E 




N 


(A.7) 


where N is total number of data points. The term (/, ) occurring in the first 

and the second gradients is evaluated by finite difference approximation for each of the 
sensitivity coefficients present in it. For example, to evaluate sensitivity coefficient 
da 


■ , the equations of motion are solved to obtain a for and 


mq 


tnq mq 
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Let the respective values of or h& a* and a , then we have 


da _ a^ -a 
2AC„ 


The reliability of the estimated parameter in terms of the Cramer-Rao bounds is 
given by the square root of the diagonal elements of ^(^ 4 ) ] ' , which is anyway 
evaluated as part of the ML algorithm. 

In the present work, for estimating the aerodynamic parameters of the artillery 
shell, the vector of unknown parameters is given by ^=[ C^,C^a>^ma’Qp J ^ 

typical case with the Z (tj and ) responses being computed by the spatial 
coordinates x,y and z respectively. 
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